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We propose a new reformulation of Yang-Mills theory in which three- and four- 
gluon self-interactions are eliminated at the price of introducing a sufficient number 
of auxiliary fields. We discuss the validity of this reformulation in the possible ap- 
plications such as dynamical gluon mass generation, color confinement and glueball 
mass calculation. We emphasize the transverse-gluon pair condensation as the ba- 
sic mechanism for dynamical mass generation. The confinement is realized as a 
consequence of a fact that the auxiliary fields become dynamical in the sense that 
they acquire the kinetic term due to quantum corrections. 



1. Introduction 

In the classical level, Yang-Mills theory ^ is a scale-invariant (actually, 
conformal invariant) field theory describing massless gauge boson. This is 
because the existence of the mass term is prohibited in the classical Yang- 
Mills theory, once we require gauge invariance. 

In the quantum level, on the other hand, this is not necessarily the case. 
To quantize the Yang-Mills theory, the procedure of gauge fixing (GF) is 
necessary. The gauge-fixed Yang-Mills theory does not prohibit the mass 
term, since the gauge invariance no longer hold. Moreover, we need to 
introduce the Faddeev-Popov (FP) ghost field in the non-Abelian gauge 
group. Consequently, the BRST invariance of the total Yang-Mills action 
is expected to hold in the quantum Yang-Mills theory instead of the gauge 
invariance. In fact, it is known ^ that a mass term of the gauge boson 
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can be introduced in the Lagrangian by adding simultaneously the (gauge- 
parameter dependent) mass term of the ghost. In this case, the BRST 
transformation should be understood as the on-shell BRST transformation. 
Here on-shell means that the Nakanishi-Lautrup field B is eliminated by 
using the equation of motion for B. However, the on-shell BRST transfor- 
mation does not have the nilpotency, unless we use the equation of motion 
for the (anti)ghost. There is an option ^ to modify the off-shell BRST 
transformation so that the mass term can be made off-shell invariant un- 
der the modified transformation by adopting an appropriate GF-I-FP term. 
However, the off-shell nilpotency is lost in the modified off-shell BRST 
transformation. No one has found the mass term which is invariant under 
the off-shell BRST transformation with nilpotency. Remember that the 
nilpotency plays the key role in the proof of the unitarity of the physical 
S-matrices.^ Therefore, the failure of nilpotency endangers maintaining the 
unitarity of the theory. In fact, it is shown that the physical space specified 
by the subsidiary condition Q'^ \phys) = using the modified BRST charge 
Q'g contains the negative norm state. These facts prevent us from pursuing 
the massive Yang-Mills theory at the tree level, i.e., Yang-Mills theory with 
a mass term in the Lagrangian level. 

In view of these, we are urged to consider the Yang-Mills theory whose 
mass is purely generated through quantum effects. Thus we are lead to 
examine the dynamical mass generation of gauge boson in Yang-Mills the- 
ory. In fact, such a possibility and its implications were investigated long 
ago in the following papers: [Eichten and Feinberg, 1974], [Smit, 1974], 
[Tye, Tomboulis and Poggio, 1975], [Fukuda and Kugo, 1978], [Gusynin 
and Miransky, 1978], [Cornwall and Papavassiliou, 1982], by using non- 
perturbative methods: Schwinger-Dyson equation, Bethe-Salpeter equa- 
tion. Effective potential. Variational technique. Cooper equation, Bogoli- 
ubov transformation, etc. However, I have no time to explain the details 
of these works. 

Here we mention recent studies of the coupled SD equations for the gluon 
and ghost propagators (in the Landau gauge) initiated by Smekal, Hauck 
and Alkofer.^'^ They revealed the surprising results: Euclidean propagators 
exhibit infrared (IR) asymptotic power-law behaviors characterized by the 
critical exponent k (1/2 < k < 1): 
• IR suppression of the gluon propagator 

DTiQ^) = ^{Qy 10 as gno. (1) 
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• IR enhancement of the ghost propagator 

Afp(Q') ^ ^(Q')-'^ T oo as Q^O. (2) 

These approximate solutions lead to 

• Existence of IR fixed point 

/3(ff2(^))tO as MiO. (3) 

• Fulfillment of a color confinement criterion due to Kugo-Ojima 

lim [g^Af p(02)]-i = 1 + y(o) = 0. (4) 

However, it should be remarked that these solutions are obtained un- 
der the specific truncation of neglecting the two-loop diagrams in the 
gluon SD equation. Nevertheless, recent lattice simulations support 
some of the above results: [Bonnet, Bowman, Leinweber and Williams, 
hep-lat/0002020], [Langfeld, Reinhardt and Gattnar, hcp-lat/0107141], 
[Alexandrou, dc Forcrand and Follana, hep-lat/0112043], [Furui and Naka- 
jima, hep-lat/0305010]. 

It is obvious that the above solutions of the SD equations do not cor- 
respond to the dynamical masss generation. However, the latest results 
of SD equations seem to suggest that k = 0.5 is allowed as a solution if 
the two-loop diagrams and hence all the diagrams are included in the SD 
equations. On the other hand, it was argued by the author ^ that the 
axiomatic consideration yields k= 1/2: 

Dt{Q'^) = const. + 0{Q'^) i const, as | q, (5) 

AMQ')-i^(Qr'/2-^^Too as Q'lO. (6) 

These solutions also satisfy (3) and (4). See Fig. 1. 

To find a non-perturbative solution with dynamical mass generation is 
an interesting question in itself. However, more important issue is to answer 
the question: What is the mechanism of generating the gluon mass? This 
is because the equation can not be solved without approximation in almost 
all cases, whereas the mechanism does not depend on the approximation 
if it is true. In fact, the above solutions of the SD equations are gauge 
dependent and not systematically improvable. 
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Figure 1. The schematic behavior of RG functions, running couphng constant and 
propagators for gluons and ghosts for Yang-Mills theory in the Lorentz-Landau 
gauge, (a) Beta function vs. Running coupling, (b) Anomalous dim. of gluon vs. 
Gluon propagator, (c) Anomalous dim. of ghost vs. FP Ghost propagator. 



In this talk, we propose the gluon pair condensation ^ as a mechanism of 
gluon mass generation and color confinement in SU{N) Yang-Mills theory. 
See Fig. 2. A purpose of this work is to construct an approximate solution 
of QCD which is blessed with both dynamical gluon mass generation and 
color confinement. 

For such Bose-Einstein (BE) condensation to take place, we need the 
net attractive force between gauge boson pair. In Yang-Mills theory, gluon 
self-interactions exist. See Fig. 3. 

In the perturbation theory, short-distance gluon-gluon potential leads to 
the net attractive force although the details depend on the total spin of 
the gluon pair. See Fig. 4. In order to see the long-distance effect realizing 
the gluon confinement, we must go beyond perturbation theory. 

For this purpose, we set up a new large N expansion based on a new 
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Figure 2. Gluon pair condensation in spin singlet and color singlet channel (center 
of mass frame). 





(a) (b) 

Figure 4. Perturbative gluon-gluon interaction to the lowest order (at tree level): 
(a) attractive interaction for s—,t—,u— channel gluon exchange, (b) attractive or 
repulsive force depending on the total spin S — Si + S2 of two-gluons . 



formulation of Yang-Mills theory. It is said that N is only one free param- 
eter in SU{N) Yang-Mills theory. In general, the large TV (or the 
expansion is 

• non-pcrturbative 

• systematically improvable 

• gauge symmetry preserving 
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• related to string theory ['t Hooft,1974] 

We also suggest a new scheme of the coupled SD equation to include all 
the diagrams. For details, see 



2. Yang-Mills theory without the gluon self-interactions 

The Euclidean SU {N^) Yang-Mills Lagrangian is given by 

^YM = -^J'l^ = \{d^A. - d^A^, + gA^, x A.f. (7) 

In what follows, we use a notation := F ■ F, {F ■ G) := F^G^ = 
2iv{FG), {F X G)^ := fABCpBQC ^ -2itv{T^[F,G]), and {F G)^ := 
dABCpBQC ^ 2ti{T^{F, G}). 

We rewrite the Yang-Mills theory into an equivalent theory in which the 
gluon self-interactions are eliminated by introducing a number of auxiliary 
fields (j), -Bj^y, G^i/, ^p^, and by choosing the parameters appropriately. 

Cym =^—^{d^A. - d^A^f + ^B^^ ■ B'^-' 

- ^S"- • *[pid^A. - d^A^) + ag{A^ x A.)] 

+ + ^<P{A, ■ An + ^G,^G^^ + ^G,^S^^ 



+ f'P-'P+f'P- {A^^An + -j-V^. ■ V^^^ + -fV,, ■ (8) 
where the three-gluon interactions have been eliminated by choosing 

p = a-\ (9) 
and the four-gluon interactions have been eliminated by choosing 



3(a^ - 1)9^ _4(l-aV 



CJG = ^^^rf^, = (3/2)(a2 - l)g\ oy = 2(1 - a^)g^. 

(10) 



Here * denotes the Hodge star (duality) operation * which is defined for 
the second-rank tensor as ^^nvpaH'"^ . 

In addition, the GF-I-FP term is given by 

£gf+fp = 5b5b -A^ + ^C- -^B-B. (11) 
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By introducing further auxiliary fields (pFP,^FP, the four-ghost self- 
interactions are eliminated as 

-Cgf+fp =^{d^A^f+C- d^d^C + (1 - 05^/. • {^^.C x C) 

- ^gA^ ■ {c X d^c) - ^xai - OgHc X C) ■ (C X C) 

(jpp (j' 

H — "I^V'FP • <^FP + ofp4>fpC -C -h a'ppifFP ■ (C ★C), 

(12) 



where 



and 



aFP = m-09\ a'pp = X^{l-i)g\ (13) 



A:=a + a', ^= ^ = ^. (14) 

Finally, we obtain the total Lagrangin ^ym — ^ym + -Cgf+fp which is 
at most quadratic in the gluon field A^: 

, '^PP a2 , '^'fP 

+ <Pfp + —^Vfp ■ fFP 

+ C- di_,d''C + aFpcpFpC ■ C + a'ppipFP -{C-kC) 

+ \a^K.^''^^A^ +A^-J^, (15) 



where we have defined JC^I^ and by 

■■=5''^ [-(1 - p'){d^5^u - d,d,) - x-'d^d,] 

■■=^t - ipd^*B^, + (1 - 05(a^C X C)^ - ig{C x d^C)^. (16) 

Here is the source for the gluon field yl^. Note that IC^^ is at most linear 

in the auxiliary field and Jj^ contains the bilinear term in the ghost and 
antighost fields. [The source term for the auxiliary field will be introduced 
later to avoid the complication.] 
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An advantage of this formulation is that we can integrate out the gluon 
field Afj, exactly: 

Seym = ^ { \b^. .B^'^ + ^^^ + ^-G^.G^^ + ^^-^+ ^V^. ■ V^'^ 

CTpp 2 '^'p p 

+ -^(Ppp + -^fFP ■ (PFP 

+ C ■ df.d'^C + aFp(t)FpC ■ C + a'ppippp -{C-kC) 

- lj^[l^^''n~^J.''} + I lnDet[/C^^'^"]. (17) 

When X = 0, ap = = a'pp and the ghost self-interaction term disappears. 
Therefore, wc do not need to introduce (ppp and (ppp in (17) which decouple 
from the theory. 

Some remarks are in order. 

1. The resulting theory has one undetermined parameter a. In the 
classical level, any value of fj reproduces the original Yang-Mills theory, 
as is trivial from the above construction. However, in the quantum level, 
this is not the case. The renormalization urges the parameter a to run 
according to the renormalization scale. 

2. For a special choice cr^ = 1, we have = = 1 and hence the 
four coefficients vanish trj. = 0, which implies that all the auxiliary fields 
(f), Gij,v, Vijtv except for B^^ decouple from the theory. Therefore, this case 
reproduces the field strength formulation of the Yang-Mills theory, which 
is identified with the deformation of the topological BF theory ^^'^^ with 
the Lagrangian: 

C-BF = \b^. ■ i?"^ - ^B^" ■ *T^,. (18) 

(The signature of p and a can bo absorbed by the redefinition of the fields 
Ap, and B^;/.) When 7^ 1, therefore, we have the scalar field ^, <p and the 
symmetric tensor fields G^^, V^j, in addition to the antisymmetric tensor 
field Bfiv. 

3. Dynamical gluon mass generation, vacuum condensate 
and glueball 

If the auxiliary (composite) field ^ develops the non-vanishing vacuum ex- 
pectation value (VEV), 

(<^) := <^o 0, (19) 
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the gluons acquire the common mass Ma'- 

Ml = #0. (20) 

We argue that the mechanism of dynamical gluon mass generation is the 
pair condensation of transverse gluons. The VEV is obtained as a non- 
trivial solution of the gap equation: 

+ \g\T(K-;},m = 0, (21) 

where is the 't Hooft coupling := g^^c and /C^J is the gluon full prop- 
agator in the condensed vacuum. We can in principle calculate the effective 
action S'£;[0] from which the gap equation is obtained as the saddle-point 
equation. The fluctuation mode of S'e[0] around c^q gives the glueball. 
For a constant 0, the effective action reduces to the effective potential, i.e., 
SeW\ — ^('/')- The 00 gives the absolute minimum of V^((/)). See Fig. 5. 



V(^) 




Figure 5. The effective potential for the auxiliary field 0. 

The ratio of the glueball mass 0^"*" to the gluon mass is given by 

Af (0++)/AfA = \/6 ^ 2.45, (22) 

in the leading order. 

Among a number of the auxiliary fields introduced above, it is the anti- 
symmetric tensor field that is most important from the viewpoint of 
confinement. This is because the anti-symmetric tensor field can couple in 
a natural way to a surface element of the two-dimensional surface whose 
boundary is the Wilson loop and this property is desirable for deriving 
the area law decay of the Wilson loop average, i.e., quark confinement. 
Although the auxiliary fields do not have the kinetic term in the begin- 
ning by definition, the kinetic terms may be generated due to quantum 
effects. Then B^i, mediates the confining force. For this scenario to be 
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successful, the existence of the condensation (jjQ, i.e., non- vanishing vacuum 
expectation value (VEV) is indispensable, since it provides with the mass 
scale which is compensated for the dimension introduced through the ki- 
netic term generated. (Note that d'^/(j)Q is dimensionless.) Without the 
condensation 0o = 0) B^j^^, can not acquire the kinetic term. 

4. Vacuum condensate of mass dimension 2 

The existence of vacuum condensate of mass dimension two, i.e., (»4^) ^ 0, 
was argued to exist recently by several groups. The composite operator 
is not gauge invariant, while V^?/' ^md JF^,^ are gauge invariant. Does such a 
vacuum condensate have any physical meaning? 

• The minimum of J?^ along the gauge orbit can have a definite physical 
meaning, as argued by [Zakharov et al.,2001] 

(5a, j d'^x^Afi^ = j d^xA^ ■ 5^^AfJ, = j d'^xA/j. ■ d^Lj{x) 

d'^xdf,A^-ui{x). (23) 



• Lattice simulations + OPE fit performed by [Boucaud et al., 
2000,2001,2002] yield the sizable value {Af,) = (l.AGeVf and the oper- 
ator product expansion of the running coupling reads 

\ -9/44 



a{p) = apert{p)y- + R[lnj^j {A^,A^)/p' . (24) 

This implies the confining potential at short distance. (It is argued how the 
instanton or monopole contributes to this condensate.) 

Thus the existence of vacuum condensate with mass dimension two 
means the existence of non-perturbative power correction even in the high- 
energy region. 

We have proposed a on-shell BRST-invariant composite operator of 
mass dimension two. It is given in the partial gauge fixing G ^ H by^^: 



O := fi-i J d'^x tvG/H 



^A^.A^' + XiCC 



(25) 



Therefore, O is invariant under the on-shell BRST (and anti-BRST) trans- 
formation when C = 5 or A = 0. In particular, in the Landau gauge 
A = 0, d^'A^ix) = 0. Hence O is off-shell BRST invariant. 
We have argued in 
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1. The on-shell BRST closed operator O {SosO = 0) is written as a 
sum of the gauge-invariant (but nonlocal) operator O' {S^O' = 0) 
and the on-shell BRST exact part: 

= 0' + 6osO", (26) 

where 

^'=""7(2^|^/'(''^^'^''^'^'^^' 

■*— gauge-invariant transverse (but nonlocal) 

^ on-shell BRST exact (27) 

and 

So^O" = J ^^-^|^[a^(fe, S)a^Hk, L) + >Ac\k)c''\k) + h.c.]. 

<— longitudinal and scalar modes cancel the ghost and antighost 

2. The VEV of O is equivalent to the VEV of the gauge invariant 
operator O' which has zero ghost number and is written in terms 
of the transverse gauge boson alone, 

{0) = {0'), 0' = Q-' fd'xJ2hA<~;\x)r. (28) 

Therefore, (O) ^ is a gauge invariant and gauge independent 
statement. 

3. The gauge invariant operator O' is nonlocal and can be non- linear 
in the non-Abclian gauge theory. The transverse gauge boson is 
both BRST-invariant and gauge invariant, i.e., <^B^phys('^) ~ 

5. Conclusion and discussion 

We have proposed an equivalent formulation of the Yang-Mills theory in 
which a sufficient number of auxiliary fields are introduced to eliminate all 
the gluon self-interactions. 

We have set up a new large A^c expansion of SU{Nc) Yang- Mills the- 
ory on a non-pcrturbativc vacTium which is stabilized by transverse gluon 
pair condensation. Then we have given the Feynman rule of the large Nc 
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expansion. This large Nc expansion is defined on a non-trivial vacuum in 
which the vacuum condensate of mass dimension two take place. A phys- 
ical meaning of the on-shell BRST invariant composite operator of mass 
dimension two is discussed. 

We determined the mass ratio of the lowest scalar glueball to the gluon: 
M{0++)/Ma = VQ = 2.45 to the leading order. This should be com- 
pared with the potential model^^ M{0++)/Ma = 2.3. For Ma = 600 ~ 
700(MeV), M(0++) = 1A7 ^ 1.72{GeV). 

To the leading order, the static potential between a pair of color charges 
is calculated. It is given by the sum of the Yukawa-type potential and the 
linear potential. The non-zero string tension is provided by the conden- 
sation. The resulting expression is very similar to that obtained in the 
MA gauge assuming the Abelian dominance. The string representation is 
also derived within the same approximation. The corresponding confining 
string is nothing but the rigid string due to Polyakov. 
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